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where  and  are  finite-difference  operators  defined  on  a  grid 
region  ft(h)  with  boundary  9ft(h)  . 

In  the  case  where  L  and  N  are  self-adjoint,  i.e. 

(1.6)  b.(x)  *  0  ,  8.j  (x)  =  0  ,  i  =  l,2,...,n 

and  a{x)  does  not  change  sign  in  ft  while  Lh  and  Nh  correspond 
to  positive  definite  matrices  one  can  attack  this  problem  via  the  appro¬ 
priate  variational  principles  (see  Weinberger  [12])- 

In  the  general  case  if  and  arise  from  a  Galerkin  approach 
and  one  has  general  L2(fi)  convergence  results  one  may  apply  the  results 
of  Vainikko  [10]  and  Osborn  [7]. 

At  first  glance  it  would  appear  that  one  could  somehow  easily  modify 
the  finite-difference  set-up  so  as  to  be  able  to  apply  the  finite-element 
theory.  This  author  was  unable  to  do  this. 

The  difficulty  centers  about  the  following  points.  Let 

T  =  L-1N  ,  Th  =  L~]Nh  . 

In  order  to  apply  the  theory  of  [  7 ]  we  require  that 
Thf  ^  Tf  Vf  e  X 

where  X  is  an  appropriate  Banach  space,  say  L2(ft)  .  But,  the  finite- 
difference  theory  is  not  designed  to  operate  on  all  of  L2(ft).  only  on 
sufficiently  smooth  functions. 

Moreover,  If  one  attempts  to  extend  the  theory  to  all  of  L2(ft)  via 
L2  projection  or  "smoothing",  one  finds  that  the  eigenvalue  problem  be- 

i 
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comes  perturbed  in  a  manner  which  is  non-trivial  when  viewed  from  all  of 
.  On  the  other  hand,  a  glance  at  the  development  of  [7]  shows 
that  all  of  the  "real  action"  takes  place  in  certain  smooth,  finite  dimen¬ 
sional  subspaces;  the  eigenspaces.  This  remark  is  the  key  to  the  develop¬ 
ment  given  in  this  report. 


'  ~^^This  work  is  motivated  by  the  desire  to  extend  the  theory  of  iterative 

^  ^  ^ 

methods  for  the  solution  of  the  algebraic  problems  associated  with  Is  de-  £•  . 

O  .  . - . -  -  v  P’ 

veloped  injjg]  to  the  non-self-adjoint  case.  See-ffij/^For  this  reason  the 

operators  nJuwi11  approximate  N  in  a  relatively  •weak1'  manner. 

In  order  to  simplify  the  presentation  we  first  consider  a  particular 


fmodel*  problem.  Following  that  development  we  discuss  the  essential  fea¬ 


tures  of  that  discussion.^ 

. 

In  section  2  we  describe  the  basic  model  problem  and  some  well-known 
facts  about  this  model  problem. 

(in  section  3  we  develop  the  convergence  theorems. 

In  section  4  we  discuss  the  ideas  in  jour  convergence  proofs  and  their 
possible  extension  to  more  general  problems. 
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2.  THE  MOOEl  PROBLEM 
Let 

(2.1 )  fi  =  {(x,y);  0  <  x,  y  < 1} 

and  let 

(2-2)  Lu  *  -(uxx+uyy)  +  a(x,y)ux  +  b(x,y)uy 

where  a,  b  e  C are  nice  functions.  Let 
(2.3)  Nu  =  a(x,y)ux  +  B(x,y)uy  +  q(x,y)u 

with  a,  6,  q  e  C°°(n)  .  For  every  f  e  L2 («)  we  let 
Tf  *  u 

be  the  solution  of 

(2.4a)  Lu  =  Nf  ,  in  , 

(2.4b)  u  *  0  ,  on  . 

Since  N  maps  L«(fi)  into  H_,(n)  and  the  resulting  u  is  in 
o  c 

H^(H)  we  see  that  T  is  a  compact  map  whose  spectrum  o(T)  consists  only 
of  eigenvalues  and  0  .  Thus  a  point  z  e  $  ,  z  t  0  ,  is  In  the  spectrum 
of  T  only  if  there  is  a  u  «  L£(fi)  which  satisfies 

(2.5)  Tu  *  zu  ,  u  $  0  . 

We  assume  that  0  is  not  an  eigenvalue  of  T  . 
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Let  P  be  a  fixed  positive  integer  and  set 

.  1 
h  “  P+1  ‘ 

Let 

(2.6a)  n(h)  =  {(x^y^)  *  (ih.jh);  l±i,j£P) 

be  the  set  of  interior  mesh  points  while 

(2.6b)  3£}(h)  i  ((x^  .y^.)  =  (ih.jh);  i  or  j  =  0  or  P  +  1 } 

is  the  set  of  boundary  mesh  points.  A  "grid  vector"  {U.  .;  0_<i,j^P  +  l} 
is  a  function  defined  on  the  total  set  of  mesh  points 

?2(h)  *  n(h)  u  3fl(h)  . 


As  usual,  we  define  the  discrete  Laplace  Operator  by:  for  i  £  i,j  £  P 

(2.7)  [a.U].  .  *  1  J - LxJ - 1-JLiJ.  +  ±ti±I - LJ - bJll 

h  i  ,j  hZ  h2 


Let 


(2.8a) 


V)  •  a(VV  •  btj  -  b(x,.y,) 


With  this  notation  we  now  define  the  finite-difference  operator  Lh  cor¬ 
responding  to  L  .  For  1  <  1,j  <  P 


(2.8) 


Remark:  While  U  Is  defined  on  the  entire  mesh  region,  the  values  A^U 
and  L^U  are  defined  only  on  the  interior  mesh  points,  n(h)  . 


VI 


w 
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We  also  define  the 
(2.9a)  CyL.j 
(2.96)  cyi^ 


usual  difference  operators 
U.+1  .-U.  . 

=  Ill  >1 — lii  ,  1  <  j  <  P,  0  <  i  <  P 
U.  -  U-  . 

=  L.JU- — U 1  ,  0  <  j  <  P,  1  <  i  <  P  . 


Let  U,V  be  grid  vectors.  We  define  an  inner  product 

P  P 

(2.10a)  (U,V)  =  h2  l  l  U.-.V 

i=1  j=l 

Observe  that  this  inner  product  involves  only  the  interior  mesh  points, 
ft(h)  .  In  the  usual  way  we  define  the  semi -norm 


(2.10b)  ||U|lj;  -  (U,U) 


We  note  that  this  semi-norm  is  a  norm  for  those  grid  vectors  U  which 
vanish  on  the  boundary,  i.e.. 


(2.11)  U  =  0  on  3fi(h)  . 


We  recall  the  following  basic  facts. 

Lemma  2.1:  Let  U  be  a  grid  vector  which  satisfies  (2.11).  Then 
(2.12)  (-AhU,U)  =  (VxU,VxU)  +  (VyU,VyU)  . 

Proof:  Apply  summation  by  parts.  See  [1],  [11]. 

Definition:  Let  f(x,y)  €  C (sT)  .  Then  Q^f  is  the  grid  vector  which 
interpolates  at  the  mesh  points.  That  is,  [Q^f]  is  the  restriction  of 


i 

* 


W"ll* 
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2  — 

Definition:  Let  f  e  C  (ft)  and  let  uJ2(f;h)  be  the  modulus  of  continuity 
of  the  second  derivatives  of  f  ,  i.e.  Let  <£>(x,y)  be  any  second  deriva¬ 
tive  of  f  ,  <p  -  fxx  or  <p  =  fyy  or  <j>  =  fxy  .  For  all  (x,y)  e  ft  and 
(x1 ,y' )  e  ft  with 

|x-x* |  <  h  ,  |y-y ' |  <  h  , 

we  have 


(2.14) 


sup  { 1 4>(x ,y )  -  $ ( x '  ,y ' )  | ;  <p  =  fxx,  f  ,  f  }  =  w2(f;h) 


Of  course,  ^(f.h)  -*■  0  as  h  0  . 

Definition:  Let  L*  be  the  formal  adjoint  of  the  operator  L  .  That  is 
(2.15a) 

Let  L*  be  the  finite-difference  operator  defined  by 


L  U  *  -(Ux/Uyy)  -  (au)x  -  (bu)y  . 


(2.15b)[LhU]ij--[AhU] 


ij 


ai+lJUi+l,.rai-ldUi-l,.i 

^i,j+lUi,j+rbid-lUi,j-l 

i  2h 

2h 

Lemma  2.2:  Let  U,V  be  grid  vectors  which  vanish  on  3ft(h),  i.e.  both  l) 
and  V  satisfy  (2.11).  Then 

(2.16)  (U,LhV)  =  (L*U,V)  . 

if  u  e  cz(n ) 


then  there  is  a  constant  K  >  0  such  that 
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(2.17a)  LhQhU  =  QhLU  +  t 

(2.17b)  L*QhU  =  QhL*U  +  t* 

where 

(2.17c)  IMI^  +  Ht*^  <  Ko)2(U,h)  . 

Proof:  The  identity  (2.16)  follows  from  an  elementary  summation  by  parts. 

The  estimate  (2.17c)  follows  from  Taylor's  theorem.  Of  course,  if  U  «  C4 

2 

then  the  truncation  error  |E|  +  |E*|  is  of  order  h  . 

Lemma  2.3:  Let  F  be  a  grid  vector  and  let 

IIHL  -  l±i,j<_P)  . 

There  are  constants  K  >  0  ,  hg  >  0  such  that,  if 

0  <  h  £  hg 


and  U  is  the 

solution  of 

(2.18a) 

L^U  =  F  in 

n(h)  , 

(2.18b) 

U  =  0  on 

3«(h)  , 

then 

(2.19) 

IIUlL  <  K||F|| 

00 

Proof:  This  well-known  result  follows  from  the  discrete  maximum  principle 

and  lemma  2.2. 

See  [11]. 
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Lemma  2.4:  Let  U  be  a  grid  vector  which  vanishes  on  3ft(h).  Then 
(2.20)  -(AhU,U)  <  |(U,LhU)|  +  Hall.  ||U||h  •  !|VxU||h 

+  iml  iiuiih  •  nyiih  • 

Proof:  Direct  Computation. 

Before  discussing  the  operators  we  discuss  the  "prolongation" 
of  grid  vectors  U  into  functions  U(x,y;h)  defined  on  ft  . 

Let  the  mesh  region  ft(h)  be  "triangulated".  To  be  specific,  imagine 
that  the  diagonal  lines 

y  =  x  -  jh  ,  j  =  -(P+1),  -  P,. . .  ,P,(P+1 ) 

have  been  drawn  over  the  grid.  Each  square  with  vertices  (x^ ,y. ), 
(xi+1,yj),  (x.+1,yj+1),  (xryj+1)  is  considered  as  the  union  of  two 

triangles. 


With  each  grid  vector  U  =  (U.,}  we  associate  a  piecewise  linear 

'  J 

function  U(x,y;h)  which  is  linear  on  each  triangle  and 


u<vyr1’1 '  uij  ■ 
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Let  P(h)  be  the  space  of  all  such  function  and  Pg(h)  be  the  sub¬ 
space  of  P(h)  which  consists  of  those  functions  which  vanish  on  3ft  . 

Throughout  the  remainder  of  this  work  we  will  make  a  complete  identi¬ 
fication  between  grid  vectors  U  =  {U^}  and  these  associated  functions 
U(x,y-,h).  This  will  avoid  cumbersome  notational  difficulties. 

Let  <  ,  >  be  the  L 2 (^)  inner  product,  i.e.  if  f,g  e  L,,(ft)  , 

then 


(2.21) 


<f,g>  = 


?(x,y)g(x,y)dxdy  . 


Lenrima  2.5:  Let  U,V  «  PQ(h)  •  Then 


<U,V>  «  (U,V)  + 


(2.22) 


2 


+  24  i^iUij(Vi+l,j-l+Vi-l,j+l) 


Furthermore 


(2.23a)  <U,V>  =  (U,V)  +  E1 


where 

(2.23b)  lE^  <  h(U,U)is(-AhV,V)h  . 

Proof:  The  formula  (2.22)  is  obtained  by  a  direct  integration.  The 
estimates  (2.23a)  and  (2.23b)  follow  from  (2.22)  and  lemma  2.1. 
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Lemma  2.6:  Let  U  €  P(h)  .  Let  Sh  be  the  triangle  with  vertices  I,J,K. 


(2.24) 


U(x,y;h)dxdy  =  -j-  (Uj+Uj-H^) 


In  this  notation  I,J  and  K  are  symbols  for  pairs  of  indices  (i,j). 
Proof:  Direct  integration. 

We  now  turn  to  a  discussion  of  the  operator  .  We  do  not  describe 
these  operators  as  specifically  as  we  described  Lh  and  L*  .  Instead, 
we  make  the  following  weak  consistency  assumption. 

o 

Assumption  N:  Let  Nh  be  defined  by:  for  1  <  i,j  <  P 


+1  ~Ui,J-l 


0  pi+1  i-ui  i  71  Pi  i 

^  ] + 


*  q(x1'yj)uij  • 

We  assume  that,  for  all  grid  vectors  U,V  which  vanish  on  3fi(h)  , 
(2.25a)  (U,NhV)  =  (U,NhV)  +  E2 


where 


(2.25b)  |  E2  j  <  Dh[(-AhU,U)-(-AhV,V)P 

+  Dh2[(-AhU,U)  +  (-AhV,V)] 

Lemma  2.7:  Let  U  €  PQ(h)  '  1”>ien  there  1S  a  constant  C.j  ,  depending  on 
ML,  ||e|L,  ||q|L,  but  not  on  h  so  that 
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(2.26a) 

(2.26b) 


o  o 


(NhU,Nhu)  <  C^-^U.U)  . 


(V,NhU)|  ^  C,  ||Vl|h(-AhU,U) 


Proof:  Direct  Computation. 


Before  proceeding  to  the  convergence  theorem  of  the  next  section,  we 
recall  some  basic  facts. 


If  U  €  Pg(h)  then 


(2.27)  (-AhU,U)  =  ||  U  ||  ^ 


H^fl) 


Moreover,  we  have  the  important 


Lemma  2.8  (Rellich):  Let  hn  -*  0  and  let  un  e  PQ(hn)  satisfy 


<-\VUn>  ■  l|U„Ho  iK 

"  H,  (B) 


Then  there  is  a  subsequence  n'  and  a  function  u  e  H1 (£2)  such  that 


Un,  —  u  in  H-j  (Q) 


U„1-’-u  in  L,(fl) 
n  i. 


Proof:  See  [2]. 
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3.  THE  CONVERGENCE  THEOREMS 


Lemma  3. 1 :  Let  [  c  |  be  a  compact  subset  of  the  resolvent  set  of  T 
Assume  that 


o  <  Z  .  . 

For  each  z  e  £  consider  the  problem 


(3.1a) 


(3.1b) 


LhU  "  zhU  =  F  *  in  * 


U  =  0  ,  on  3Q(h) 


There  are  constants  h-j  >  0  ,  K  >  0  with  K  independent  of  h  ,  such 
that:  For  0  <  h  £  h-| ,  problem  3.1  has  a  unique  solution 


U  ,  and 


(3.2) 


|U||h  <  R||F||h  . 


Moreover,  for  every  u  e  C^(fi)  which  vanishes  on  3fi,  let  V  be  the 
grid  vector  which  satisfies 


(3.3a) 


(3.3b) 


LhV  '  zNhV  =  VLu'zNu^  *  in 


V  =  0  ,  on  3fi(h)  . 


<  V-u,V-u >  -*■  0  as  h  -*•  0 


That  is,  V(x,y;h)  converges  to  u  in  Lj^) 
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Proof:  It  is  only  necessary  to  verify  (3.2).  Suppose  (3.2)  is  false. 
Then  there  are  sequences  zn  e  l  ,  hp  ■*  0  and  grid  vectors  Fn  such 


that 

(3.4a) 


l|Fnllh  *0  «  hn-0  , 


while  the  corresponding  solutions  [(3.1a),  (3.1b],  Vn,  satisfy 


Applying  lemma  2.4  we  have 

(3.6)  -uh  V„,Vn)  <  |(Vn,fn)|  +  l|a||J|VxVnl|  *  UbllJIyA 

n  nr 

Using  Assumption  N  and  Lemma  2.7  we  have 

Vn)|  <  C,(-4h  yv/  *  2»n(-V«-V») 

n  n 

Therefore 

(3.7)  -U„  V„.*„)  i  II 'A  ♦  *2  ♦  *„ -V 


(3.4b) 

n»A  * '  • 

n 

Then 

i 

♦ 

(3.5a) 

Lh  \  -  rVn  *  • 

n  n 

in  fl(hn)  , 

(3.5b) 

v„  =  o  , 

1 

on  3f2(hn)  . 

(3.5c) 

ll*„llh  ■ '  • 

1 


where 


*2  '  10<ll«ll£  *  Itblll  ♦  Cf) 

C1 

ci  =  max  {]TT  *  z  e  D  • 

5  -  max  {  -njr  ,  z  e  £}  . 


Thus,  when 


Dhn 


we  obtain 


(3.8) 


(-Ah  Vn,Vn)  <  2(1  +K2) 


Let  zn,  be  a  subsequence  of  the  zp  such  that 


v  - 5  €  i  • 


Applying  Lemma  2.8  we  see  that  there  is  a  sub-subsequence  n"  and  a  func- 
o 

tion  v  e  H1  such  that 

o 

(i )  Vn„  -  v  in  (Q) 


(11)  Vn"  -  v  in  L2(n)  . 


For  simplicity  of  notation  we  denote  {hn„}  by  (hn>  . 


Let  ip  e  Cg(n)  and  let 


(3.9) 


0n  •  Q„  0  . 
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From  (3.5a),  (3.5b)  we  obtain 


<VLh  V  -r'<VV„>  •  (Vfn> 

n  n 


Thus,  using  Assumption  N  we  see  that 


where 


<L!LvV  -f  <VV„>  *  <VV 

n  n 


|gnllh  -  0  as  n  -  " 

n 


On  letting  n  «•  lemma  2.2  and  lemma  2.5  yield 
(3.10)  <  L*4>,v  >  -  t<«J>,Nv>  *  0  . 

S 

Thus,  v  is  a  weak  solution  of 


(3.11a) 

(3.11b) 


Lv  -  rNv  *  0  ,  in  Q  , 


v  =  0  ,  on  3i2 


But,  then,  v  is  a  genuine  solution  of  (3.11a),  (3.11b)  and  ;  t  J  is 
an  eigenvalue  of  T  (see  [2]).  With  this  contradiction  we  have  proven 
(3.2). 

The  remainder  of  the  lemma  follows  from  similar  arguments. 

Let  u  t  0  be  an  eigenvalue  of  T  .  Let 

S  =  ker  (u-T)3  , 

i.e. ,  the  kernel  or  nullspace  of  the  linear  map  (u-T)3  where  a 
is  the  ascent  or  rank  of  y  (see  [5]).  S  is  the  range  of  the  linear 


’  i  k  mi 


17 


projection  E  given  by 

(3-12)  E  -  27T  \  (z_T)'1dz 

with  r  any  circle  in  the  resolvent  set  of  T  enclosing  u  in  its  in¬ 
terior,  but  no  other  point  of  the  spectrum.  The  linear  projector  E  can 
also  be  written  in  the  form 

(3.13a)  Ey  =  f  <y,4>?><b. 

i=l  1  1 

where  {<)>.,  i=l,2,...a)  is  any  basis  for  S  and  i=l,2...a)  is  the 

corresponding  dual  basis  in 

(3.13b)  S*  *  ker  (u-T*)a  . 

That  is  «  S*  and  =  6.^  .  Here  T*  is  determined  by 

(3.13c)  <Tu,v>  =  <u,T*v>  . 

That  is, 

(3.13d)  T*  *  N*(L*)_1  . 

Let  l  *  r  .  Applying  Lemma  3.1  we  see  that 

(3.14)  U-Thr’  *I(Lh-JNh>',Lh 

is  well-defined  on  r  for  0  <  h  <_  h^  .  Let 

<3-,5>  Eh  *  5J7 1  • 

1 
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Lemma  3.2:  Let  u  t  0  be  an  eigenvalue  of  T  .  There  is  an  >  0  such 
that  if  0  <  h  <  h2  there  is  an  eigenvalue  uh  of  Th  which  lies  in  the 
interior  of  r  . 

Proof:  If  the  lemma  is  false,  then 
EhE°. 

Let  <{>  £  S  be  an  eigenfunction  of  T  with  eigenvalue  u  •  It  is  not  dif- 

ficult  to  prove  that  4>  e  C  (ft)  .  However,  in  order  to  see  how  the  general 

o 

argument  goes,  we  replace  by  a  smooth  H2  n  ^  approximant,  which  we 
again  call  <j>  .  Let 

(3.16a)  *hsQh*.  *h*Eh0h 

and  let 

(3.16b)  ¥h  =  QhE$  , 

we  observe  that 


fh(x,y;h)  -  E  f  0  . 

We  calculate 

’h  -  *h  ■  ST  j  f  L‘-  ViV'Vi,’  • 

r 

Since  <$>  is  smooth  we  may  apply  lemma  2.2.  Thus 

Lh‘h  ’  v*  * 1 
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where  IMi^  -*■  0  .  In  fact,  for  each  fixed  z  e  T  let  v(x,y;z)  be  the 
solution  of 

(3.17a)  (L-jN)v  *  L<}>  ,  v  =  0  on  3ft  , 

while  V  is  the  solution  of 

(3.17b)  (Lh'  zVV  =  Vh*  ’  V  =  0  on  3ft(h). 

From  Lemma  3.1  we  see  that 


HVIIh  *  (-V-V)11  <  (C||U||„  . 

Thus,  using  Lemma  2.5  and  Holder's  inequality  we  see  that  there  is  a  constant 
K1  such  that 


(3.18) 


|V||L2(Q)  -  K1  1  l,V,lL1  («)  -  K1 


Thus 


(3.19) 


u(x,y;h)  -  V(x,y;h)  * 


~  (v(x,y;z)  -  V(x,y;z))dz 


Applying  Lenina  3.1  we  see  that 


2ir 


u-V||  2  —  2tF  0d6  l|v(*,-;z)-V(.,.;z)||  2 

LZ(ft)  Z7r  L  LZ(ft) 


However,  according  to  Lemma  3.1  the  integrand  goes  to  0  for  every  z  £  r 
and  the  estimates  (3.18)  allow  us  to  apply  the  dominated  convergence  theorem. 
Thus 

V(x;y;h)  -  E*  t  0  . 


I 


Corollary:  Let  denote  the  range  of  Eh  .  Then,  for  0  <  h  £  h.. 


we  have 


dim  Sh  _>  dim  S  =  a  . 


Proof:  We  need  only  apply  the  above  argument  to  each  <{>..  in  the  basis  of 
S  ,  smoothed  if  necessary.  The  resulting  grid  vectors  E^Q^.  will  be 
linearly  independent  if  h  is  small  enough. 


Lemma  3.3:  dim  S  =  dim  Sh  . 


Proof:  Let  y^  f  0  and  be  an  eigen  pair  of  T^  .  That  is. 


LuVk  =  —  N  hV. 
h  h  uh  h  h 


Suppose  y^  -*•  y  .  Then,  it  follows  from  Lemma  3.1  that  y  is  an  eigenvalue 


of  T 


In  fact,  let  V,  satisfy 


<vV  \  -  \  * 0 


(vV  Vh  .  o  . 

A  modification  of  the  argument  of  Lemma  3.1  shows  that,  after  selection  of 
a  subsequence  if  necessary,  there  are  functions  v,w  e  L2(n)  such  that 


Vh  ■*  v  ,  -*■  w  +  0 


(y-T)r"^v  =  w  f  0 


(y-T)  v  *  (y-T)w  =  0 


Thus,  a  basis  for  Sh  converges  to  a  set  of  linearly  independent  func¬ 
tions  in  S  .  Therefore 


Urn  sup  dim  Sh  <_  dim  S  =  a  . 

Thus,  the  lemma  follows  from  the  corollary  to  Lemma  3.2. 

We  collect  our  results  in  the  following 

Theorem  3.1:  Let  y  f  0  be  an  eigenvalue  of  T  with  algebraic  multi¬ 
plicity  a  .  Let  r  be  a  circle  about  y  which  does  not  include  any  other 

point  of  the  spectrum  of  T  either  in  its  interior  or  on  its  boundary.  Then 

there  is  ah  h1  >  0  such  that  if  0  <  h  £  h^ 

(i)  there  are  no  eigenvalues  uh  of  T^  on  r 

(ii)  there  are  eigenvalues  u1>h,  y2>h,...yr  h  of  Th  which  lie 
inside  r  . 

(iii)  Let  m^,m2,...,mr  be  the  algebraic  multiplicities  of  y^  h, 
u2  h,...yr  h  respectively. 


Then 


Then  is  also  one  to  one  and  onto  from  to  S  . 


Once  one  has  this  result  one  can  follow  Kreiss  [6],  Osborn  [7]  or 
de  boor  and  Swartz  [3]  to  obtain  Theorem  3.2:  Let 


1 


u  =  —  y  m.u.  t 
a  J  J.h 


Then  there  is  a  constant  K  such  that 


|v-u(  1  K  sup  (T-Th) |  . 

Specifically,  we  have 

Let  ,  i  =  l,2,...a  be  a  basis  for  S  and  (/>T,1  =1,2... a  be  a  basis 
for  S*  .  Then  there  is  a  constant  K  such  that 

ju-u!<K  l  <(T-Th)<i>a,4,*> 
a=l 

Proof:  See  the  development  in  [3]. 


4.  THE  GENERAL  IDEA 

While  the  discussion  in  section  3  is  limited  to  the  model  problem 
described  in  section  2,  one  can  see  the  essential  points  and  describe  the 
basic  ideas  which  will  enable  one  to  discuss  more  general  problems. 

The  basic  tool  is  the  relationship  between  the  linear  projectors  E 
given  by  (3.12)  and  Eh  given  by  (3.15). 

A  look  at  the  proof  of  Lemma  3.2  shows  that  one  need  only  apply  Eh 
to  "smooth"  elements  4>h  .  However,  even  when  we  apply  to  a  smooth 
,  we  must  have  results  such  as  those  given  in  Lemma  3.1  . 

Thus  it  is  worthwhile  to  study  the  argument  in  Lemma  3.1.  Almost  any 
discussion  of  the  operator  L^  yields  estimates  such  as  those  given  in 
Lemma  2.2  (consistency)  and  Lemma  2.3  (convergence  in  the  smooth  case). 
However,  the  use  of  E  and  Eh  requires  estimates  on  the  solutions 
of  the  equations 

zLhV  -  Nhv  =  fh  ,  z  e  r  . 

We  have  essentially  used  estimates  on  L^  to  obtain  these  estimates.  The 
basic  facts  are 

FACT  1:  If  {V,.;  0<h<hn}  is  the  set  of  all  "functions"  in  P(h)  such 
-  n  —  u 

that 

l|LhvhllhiC,  0<hih0, 
then  the  (Vh>  are  sequentially  compact. 
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FACT  2:  If  fh>  Vh  e  P(h)  and 

<4-la)  2LhVh'  Vh  =  fh 

(4.1b)  fh  -  f  e  C2 

(4.1c)  Vh  -  V 

then,  in  an  appropriate  sense 
(4.2)  zLV  -  NV  =  f  . 

In  the  proof  of  Lemma  3.1  we  first  use  these  facts  to  establish  (3.2). 
Then  we  use  these  facts  and  (3.2)  to  prove  the  convergence  of  the  solutions 
of  (3.3a),  (3.3b). 

Thus,  in  the  more  general  case,  i.e.  general  L,  N  and  fi,  one  would 
expect  to  prove  analogous  spectral  convergence  theorems  if  one  can  establish 
these  two  facts.  And,  our  discussion  shows  that  this  is  the  case.  The  proof 
of  these  facts  are,  in  general,  not  too  difficult.  We  have  chosen  to  use 
"energy  estimates"  to  establish  FACT  1  and  Weyl's  lemma  (weak  solutions  are 
genuine  solution)  to  establish  FACT  2. 

In  any  particular  instance  one  has  a  variety  of  possible  arguments  which 
lead  to  such  conclusions.  For  example,  in  the  case  discussed  in  this  report, 
i.e.  the  model  problem,  one  could  use  the  estimates  of  A.  Brandt  [4]  to  es¬ 
tablish  FACT  1.  In  fact,  in  this  case,  one  can  use  the  Banach  space  C(fl) 
rather  than  Lg^)  and  rely  completely  on  the  estimates  of  [4].  We  choose 
to  use  "energy  estimates"  because  we  feel  that  they  are  (i)  weaker  and  (ii) 
(generally)  easier  to  obtain. 
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